Modulational instability in nonlocal Kerr-type media with random parameters 
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Modulational instability of continuous waves in nonlocal focusing and defocusing Kerr media with 
stochastically varying diffraction (dispersion) and nonlinearity coefficients is studied both analyt- 
ically and numerically, ft is shown that nonlocality with the sign-definite Fourier images of the 
medium response functions suppresses considerably the growth rate peak and bandwidth of insta- 
bility caused by stochasticity. Contrary, nonlocality can enhance modulational instability growth 
for a response function with negative-sign bands. 



PACS numbers: 42.25.Dd, 42.70.-a, 42.65.Jx 
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I. INTRODUCTION 

Modulational instability (MI) in nonlinear media is a 
destabilization mechanism which produces a self-induced 
breakup of an initially continuous wave into localized 
(solitary wave) structures. This phenomenon was pre- 
dicted in plasma [U, 0] , nonlinear optics^jHjfluids Q 
and atomic Bose-Einstein condensates [a, H Q. MI of 
continuous waves can be used to generate ultra-high 
repetition-rate trains of soliton-like pulses [1, [13, lU'l • It is 
common knowledge that MI is absent in the defocusing 
Kerr medium and presents as the long-wave instab ility 
with a finite bandwidth in the focusing Kerr medium [ij] • 

The above results were obtained for media with de- 
terministic parameters. Contrary, in realistic media the 
characteristic parameters are not constants, as a rule, 
but fluctuate randomly around their mean values. It was 
shown in the setting of nonlinear optics that stochastic 
inhomogeneities in a Kerr-type medium extend the do- 
main of MI of continuous waves, as compared with deter- 
ministic systems, over the whole spectrum of modulation 
wavenumbers, even for the defocusing regime p^. [Tsj . 
A comprehensive review of MI of electromagnetic waves 
in inhomogeneous and in discrete media is given in 
Ref. ''16]. 

Another important aspect of a class of realistic nonlin- 
ear media is concerned with their nonlocality. Nonlocal- 
ity is typically a result of underlying transport processes 
such as heat conduction in thermal nonlinear media [l3l , 
diffusion of atoms in a gas I18| . long-range electrostatic 
interaction in liquid crvstals |l9|. charge carrier transfer 
in photorefractive crystals j2Cll. \2l\. and many- body in- 
teraction in Bose-Einstein condensates [22]. Nonlocality 
can prevent the collapse of self-focused beams 23, 2^] and 
dramatically alter interaction between dark solitons [1^ . 
MI in deterministic nonlocal Kerr-type media was stud- 
ied in Refs. [26, 27], and it was shown that nonlocality 
does not produce MI in the defocusing case for small and 
moderate values of the product "modulation amplitude 



X nonlocality parameter" . 

In the present paper we unite the two above lines of 
study of nonlinear media and analyze MI in nonlocal me- 
dia with stochastic parameters. Since nonlocality spreads 
out localized excitations, it is reasonable to expect a par- 
tial suppression of the stochasticity-induced MI gain. In- 
deed, we demonstrate that the aforementioned situation 
with MI in local stochastic media with the sign-definite 
Fourier images of the response functions changes drasti- 
cally, if nonlocality is taken into account. Namely, both 
the growth rate peaks and bandwidths of instability are 
considerably decreased. On the other hand, there can be 
an "anomalous" behavior of nonlocality when the Fourier 
image of the response function of a nonlocal medium al- 
lows for sign-negative bands. In this case the MI gain of 
a nonlocal medium can exceed that of a local stochastic 
medium for some values of the modulation wavenumber. 
We adopt the nonlocal nonlinear Schrodinger equation 
with random coefficients as a model to reveal peculiar- 
ities of MI of continuous waves. The results obtained 
are illustrated by the white noise model for parameter 
fluctuations and by response functions of several types. 



II. MODEL 

The propagation of an optical beam along the z axis in 
a nonlocal medium with random parameters is governed 
by the nonlinear Schrodinger equation 
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+ g{z)u / dx' R{x — x' )\u\'^ {x , z) = . 

(2.1) 

Here x is the transverse coordinate, u{x, z) is the complex 
envelope amplitude and we use the standard dimension- 
less variables. The group velocity dispersion (or diffrac- 
tion) coefficient d{z) and nonlinearity coefficient g{z) are 
considered as stochastic functions which fluctuate around 
their mean values dp (do > 0) and go (go ^ 0): 

d{z)=da{l + md{z)), g{z) = go{l + mg{z)). (2.2) 

Here and mg are independent zero-mean random pro- 
cesses of the Gaussian white-noise type, 

(md) = (mg) = 0, (■md{z)md{z')) = 2(t^(5(z - z'), 
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{m,{z)m,{z'))=2al5{z-z')), 

and the angle brackets stand for the expectation with 
respect to the distribution of the processes md{z) and 
mg{z). The integral in equation (12. ip represents the field- 
intensity dependent change of the refractive index char- 
acterized by the normalized symmetric response function 
R{x)^ J_^dxR{x) — 1. The delta-function response 
function R{x) = 6{x) corresponds to the local limit of 
the model. We will discriminate between the focusing 
(50 > 0) and defocusing {go < 0) media. 

Eq. (j2.1|) possesses the homogeneous plane wave solu- 
tion 



uo = A exp 



lA' I dz'g{z') 




(2.3) 



where A is a real amplitude. Now we perform the linear 
stability analysis of the solution (|2.3[) . Assume that 



If we were deal with the deterministic system with the 
parameters d n an d go, Eq. (|2.6p would be the main object 
to study MI [2^. However, MI induced by the random 
fluctuations is not captured by the analysis of the first 
moments (p) and (a) |15| , and it is necessary to compute 
the modulational intensity growth given by the higher- 
order moments. 

III. THE SECOND-ORDER MOMENT MI GAIN 

We consider the second moments (p^), (per) and (u^) 
as constituents of the column vector 



(3.1) 



The moment (per) is added to close the equations for 
the second-order moments. Then we should calculate z- 
evolution of the vector X^'^\ Its first component gives 



m(x, z) = {A + u(x, z)) exp 



lA^ / dz'g{z') 



(2.4) 



is a perturbed solution of Eq. (I2.1|l with v{x, z) being a 
small complex modulation. Substituting Eq. (|2.4p into 
Eq. (|2.1[) and linearizing about the plane wave (|2.3p . we 
get a linear equation for v{x, z): 

iv^ + ^d{z)vxx + 2g{z)A'^ J dx' R{x - x')Re v{x' , z) = 0. 

(2.5) 

After decomposing v into real and imaginary parts, v = 
r(x, z) -\- is{x, z), and performing the Fourier transforms 



p{k,z) = 



1 f°° 

— J dxr{x,z)e^^, 



-r{P^) = 2(p^p) = doP{pa) + doP {md{z) pa) , 
dz 



in accordance with Eqs. (|2.2p and (|2.6p . For decoupling 
of the mean (md(z)pa) we apply the Furutsu-Novikov 
formula [28I. [29I 



{md{z)pa) = / dyajBiz - y) 



Sm, 



-pa 



(3.2) 



Here B{z — y) = 6{z — y) for the white-noise Gaussian 
random process, while the functional derivative (S/Srud) 
is calculated from Eq. (|2.6p . Indeed, writing p{z) as the 
integral p{z) = {l/2)k'^ dyd{y)a{y) (and the similar 
integral for a{z)) and accounting for the explicit repre- 
sentation (|2.2p of d{z) in terms of gives 



a(k, z) 



2ti 



dx s{x,z)e'-'"' , 



SjP'^) ^ JP_^ ^ ^ -dok'^ia'^ 

Snid Snid Snid 2 



R{k) = ^ / da;i?(x)e^'=^. 



Eq. (|2.5p is converted to a system of linear equations for 
p and a: 



d_f p 
dz\ cr 







-^d{z)P + 2g{z)A^R 



\d{z)k'^ 



(2.6) 



Therefore, {md{z)pa) = {l/2)a^dok^{{a^) - (p^)) and 
finally 



±y) ^ doPipa) + lajdlk^i{a-') - ip^)). 



Just in the same way we can calculate z-derivatives of 
the other components of the vector X'^^K As a result, we 
obtain the evolution equation {d/dz)X^'^'> = Mf^)^^^) 
with the 3x3 matrix M^^^ of the form 
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- 1 



2[2goA^R--dok^ 



dok^ 



(3.3) 



Eigenvalues of M^^^ with positive real parts lead to in- 
stabilities, and the largest positive value determines the 
MI gain G2 (k) . The eigenvalues Xj are easily found from 
Eq. (|3.3p but they are too cumbersome to be reproduced 
here explicitly. Below we separately analyze the cases of 
the defocusing {go < 0) and focusing {go > 0) nonhn- 
earities. Following [2^, we will use for illustration the 
Gaussian response function 



Rg{x) 



1 I X 

:exp ( - — 



and the exponential one 



i?G(fc) = exp ( --a^fc^ 
(3.4) 



Re{x) 



2a 



exp 



Re{k) 



1 



(3.5) 



as examples of the response functions with the sign- 
definite Fourier images, as well as the rectangular re- 
sponse function 



Rr{x) 



— for Ixl < a, 
2a II-' 



for Ixl > a. 



Rr{k) 



sm{ak) 
ak 



(3.6) 

whose Fourier transform has negative-sign bands. Here a 
is the nonlocality parameter, a — *■ means R{x) — > 5{x) 
and R{k) 1. 




FIG. 1: Defocusing media. Plots of the MI gain G2{k) for 
local stochastic medium (solid line), nonlocal stochastic media 
with the Gaussian (dash-dotted line) and exponential (dashed 
line) response functions. Here do = 2, |(?o|^^ = 1, 
0.1. Upper panel: a — 1; lower panel: a = \/2. 



A. Defocusing nonlinearity 

For the defocusing nonlinearity go < we obtain one 
real eigenvalue Ai and two complex conjugate ones A2 
and A3. Numerical analysis shows that Ai is positive for 
all fc2, while A2 and A3 have negative real parts for the 
Gaussian and exponential response functions. Let us re- 
mind that there is no MI for go < for local deterministic 
Kerr media, while randomness of the coefficients d{z) and 
g{z) completely destroys stability of the continuous wave 
solution. This situation considerably changes for nonlo- 
cal media. Indeed, Fig. [1] clearly shows that nonlocality 
with the sign-definite response functions suppresses both 
the growth rate peak of G2{k) = Ai and MI bandwidth, 
the latter being practically finite. When the nonlocality 
parameter a grows, the suppression effect becomes more 
pronounced. Somewhat different situation takes place 



for the rectangular response function (|3.6p . For sufR- 
ciently high nonlocality, MI gain maximum for a given 
wavenumber k can exceed the corresponding value of G2 
for a local random medium (Fig. [2|). Besides, the MI 
bandwidth becomes strictly finite in this limit. 



B. Focusing nonlinearity 

In the case of the focusing nonlinearity {go > 0) a local 
deterministic medium produces the long-wave instability 
with a finite bandwidth. Stochasticity of medium param- 
eters extends the bandwidth to the whole spectrum of 
modulation wavenumbers. Calculation of eigenvalues of 
the matrix M2 ()3.3|) with go > demonstrates that non- 
locality suppresses the MI gain and bandwidth for media 
with both sigh-definite (Fig. [3]) and sign-indefinite (Fig. 
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FIG. 2: Defocusing media with the rectangular response func- 
tion. Plots of the MI gain G2(k) for local stochastic medium 
(solid line), nonlocal stochastic media with a = 2 (dashed 
line), a = 6 (dash-dotted line), and a = IQ (dotted line). 
Other parameters are the same as in Fig. [1] 



FIG. 4: Focusing media. Plots of the MI gain G2{k) for a local 
deterministic medium (solid line), local stochastic medium 
(dotted line), nonlocal stochastic media with the rectangular 
response function: a = 2 (dash-dotted line), a = 8 (dashed 
line). Other parameters are the same as in Fig. |3] 
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FIG. 3: Focusing media. Plots of the MI gain G2{k) for a lo- 
cal deterministic medium (solid line), local stochastic medium 
(dotted line), nonlocal stochastic media with the Gaussian 
(dash-dotted line) and exponential (dashed line) response 
functions. Here do = 2, goA^ = 1, crj = 0.1, = 0.2. 
Upper panel: a — 1; lower panel: a = 2.5. 



S]) response functions. Notice that stronger nonlocality 
is needed for focusing media to achieve a reduction of 
the MI gain, as compared with defocusing ones. Besides, 
maximum positions of the MI gains shift toward smaller 
wavenumbers k under nonlocaUty growth, producing fi- 
nite bandwidth. 



IV. HIGHER-ORDER MOMENTS 

The second-order moments (j3.1[) do not provide an 
analysis of the MI gain in stochastic media with sufficient 
detail. In particular, it is important to see fluctuations 
of the exponential growth of the modulation amplitude. 
More deep insight into the problem demands to account 
for higher-order moments 

= |< j = 0,...,2n. (4.1) 

In this section we study the interplay of nonlocality and 
exponential growth of the higher moments X^^"^ in virtue 
of stochasticity. As before, applying the Furutsu-Novikov 
formula p.2p . we obtain a matrix M^^"' in the form 

M(2") = dofc2^(2") + (2^0^^^? (4.2) 

2 

Non-zero entries of the matrices A'-^"\ B^^"-\ C(2") ^nd 

Then the maximal real part of roots of the characteristic 
polynomial det |M(2n) _ ^/| gjye nG2n{k)- Since all 
the matrix elements of M'^2"-) are real and the character- 
istic polynomial is of the odd degree, at least one of the 
eigenvalues of M(2") is real and the others are mutually 
complex conjugate. In what follows we will consider the 
4-th and 6-th moments. 



5 



A. Defocusing nonlinearity 

In Fig. [5] we show the results of calculating MI gains 
G2 , G4 and Gq for both the exponential and Gaussian re- 
sponse functions and compare them with the same curves 
for local stochastic media obtained in [15]. It is seen 
that nonlocality suppresses the higher-order moments as 
well. Notice that in defocusing media positions of MI 
gain maxima for moments of different orders coincide (isj 
(they are deterministic rather than random). Nonlocality 
does not disturb this property. Fig. [6] demonstrate simi- 
lar curves for the rectangular response function for differ- 
ent values of the nonlocality parameter a. It is seen that 
for sufhciently high a the medium demonstrates prac- 
tically coinciding distributions of higher-moment growth 
rates, their maxima being shifted to shorter wavelengths. 
Evidently, higher-order moments for the rectangular re- 
sponse function manifest the same "anomalous" enhance- 
ment of the growth rate in a narrow region of modulation 
wavenumbers, as compared with the local stochastic case. 
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FIG. 5; Defocusing media. Plots of the MI gains Ge (solid 
line), G4 (dashed line), and G2 (dash-dotted line) for a lo- 
cal stochastic medium (upper three curves) and for nonlocal 
stochastic media (lower three curves). Here do = 2, — 1, 
|go|j4^ = 1, (tJ = (Tg = 0.1. Upper panel: exponential re- 
sponse function; lower panel: Gaussian response function. 



B. Focusing nonlinearity 

For the focusing media the higher-order MI gains 
demonstrate much the same behavior as for the defo- 
cusing ones for both sign-definite and sign-indefinite re- 
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FIG. 6: Defocusing nonlocal stochastic media with the rect- 
angular response function. Plots of the MI gains Ge (solid 
Une), Gi (dashed line), and G2 (dash-dotted line) for do ~ 2, 
Igol^^ = 1, a J = ag = 0.1. Upper panel: a — 2; middle panel: 
a = 6; lower panel: a = 8. 



sponse functions. Fig. [7] shows the MI gains for local 
and nonlocal stochastic media for the Gaussian response 
function. Curves for the exponential and rectangular re- 
sponse functions are qualitatively the same. With in- 
creasing the nonlocality parameter a, curves for MI gains 
of different orders become closer one the other, so high 
nonlocality smoothes fluctuations of the modulation am- 
plitude growth. 
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FIG. 7: Focusing media. Plots of the MI gains Ge (solid 

line), G4 (dashed line), and G2 (dash-dotted line) for a lo- 
cal stochastic medium (upper three curves) and for nonlo- 
cal stochastic medium with the Gaussian response function 
(lower three curves). Here do = 2, a = 2, goA^ = 1, 
al = a^„= 0.1. 



CONCLUSION 



nonlocal nonlinear Kcrr-type medium with random pa- 
rameters. For the case of the white-noise model of pa- 
rameter fluctuations, we derived the equations which gov- 
ern the dependence of the MI gain on the modulation 
wavenumber. As was expected from physical motiva- 
tions, nonlocality causes c;onsidcrable suppression of the 
stochasticity-induced MI growth rate for media with the 
sign-definite Fourier images of the response functions. At 
the same time, nonlocal media with the sign-indefinite 
Fourier images of the response functions can display a 
somewhat diff'erent behavior leading to an increase, as 
compared with local media, of the MI gain for some do- 
mains of modulation wavenumbers. 



Acknowledgments 



Within the limits of the linear stability analysis, we 
have investigated the MI of a homogeneous wave in a 



The authors are very grateful to F. AbduUaev and J. 
Garnier for constructive comments. 



[1] G.A. Askar'yan, Zh. Eksp. Teor. Fiz. 42, 1576 (1962) [16 

[Sov. Phys. JETP 15, 1088 (1962)]. 
[2] T. Taniuti and H. Washini, Phys. Rev. Lett. 21, 209 [17 

(1968). 

[3] V.I. Bespalov and V.I. Talanov, Pis'ma Zh. Eksp. Teor. [18 

Fiz. 3, 471 (1966) [JETP Lett. 3, 307 (1966)]. [19 
[4] L.A. Ostrovskii, Zh. Eksp. Teor. Fiz. 51, 1189 (1966) 

[Sov. Phys. JETP 24, 797 (1967)]. [20 
[5] T.B. Benjamin and J.E. Fair, J. Fluid Mech. 27, 417 

(1967). [21 
[6] B. Wu and Q. Niu, Phys. Rev. A 64, 061603 (2001). 
[7] V.V. Konotop and M. Salerno, Phys. Rev. A 65, 

021602(R) (2001). [22 
[8] L. Salasnich, A. Parola, and L. Reatto, Phys. Rev. Lett. 

91, 080405 (2003). [23 
[9] K. Tai, A. Hasegawa, and A. Tomita, Phys. Rev. Lett. 

56, 135 (1986). [24 
[10] G. Millot, E. Seve, S. Wabnitz, and M. Haelterman, J. 

Opt. Soc. Am. B 15, 1266 (1998). [25 
[11] K.E. Steeker, G.B. Partridge, A.G. Truscott, and 

R.J. Hulet, Nature 417, 150 (2002). [26 
[12] Yu.S. Kivshar and G.P. Agraval, Optical Solitons: From 

Fibers to Photonic Crystals, Academic Press, San Diego, [27 

2003. 

[13] F.Kh. Abdullaev, S.A. Darmanyan, S. Bishoff, and 

H.P. S0rensen, J. Opt. Soc. Am. B 14, 27 (1997). [28 
[14] M. Karlsson., J. Opt. Soc. Am. B 15, 2269 (1998). [29 
[15] J. Garnier and F.Kh. Abdullaev, Physica D 145, 65 
(2000). 



F. Kh. Abdullaev, S.A. Darmanyan, and J. Garnier, 

Progr. Opt. 44, 303 (2002). 

A. Dreischuh, G.G. Paulus, F. Zacher, F. Grasbon, and 
H. Walther, Phys. Rev. E 60, 6111 (1999). 

D. Suter and T. Blasberg, Phys. Rev. A 48, 4583 (1993). 
C. Conti, M. Peccianti, and G. Assanto, Phys. Rev. Lett. 
91. 073901 (1993); 92,113902 (2004). 

M. Segev, B. Crosignani, A. Yariv, and B. Fischer, Phys. 
Rev. Lett. 68, 923 (1992). 

G. C. Duree, J.L. Shultz, G.J. Salamo, M. Segev, 
A. Yariv, B. Crosignani, P. Di Porto, E.J. Sharp, and 
R.R. Neurgaonkar, Phys. Rev. Lett. 71, 533 (1993). 
V.M. Perez-Garcia, V.V. Konotop, and J.J. Garciar 
RipoU, Phys. Rev. E 62, 4300 (2000). 

S.K. Turitsyn, Teor. Mat. Fiz. 64, 226 (1985) [Theor. 
Math. Phys. 64 797 (1985)]. 

O. Bang, W. Krolikowski, J. Wyller, and J.J. Rasmusscn, 
Phys. Rev. E 66, 046619 (2002). 

A. Dreischuh, D. Neshev, D.E. Peterson, O. Bang, and 
W. Krolikowski, Phys. Rev. Lett. 96, 043901 (2006). 
W. Krolikowski, O. Bang, J.J. Rasmussen, and J. Wyller, 
Phys. Rev. E 64, 016612 (2001). 

W. Krolikowski, O. Bang, N.J. Nikolov, D.N. Neshev, 
J.J. Rasmussen, J. Wyller, and D. Edmundson, J. Opt. 
B: Quantum Semiclass. Opt. 6, S288 (2004). 
K. Furutsu, J. Res. NBS, D-67, 303 (1983). 

E. A. Novikov, Sov. Phys. JETP 20, 1290 (1964). 



